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We demonstrate for the first time that three sterile neutrinos alone can simultaneously explain 
neutrino oscillations, the observed dark matter and the baryon asymmetry of the Universe without 
new physics above the Fermi scale. The key new point of our analysis is leptogenesis after sphaleron 
freeze-out, which leads to resonant dark matter production, evading thus the constraints on sterile 
neutrino dark matter from structure formation and x-ray searches. We identify the range of sterile 
neutrino properties that is consistent with all known constraints. We find a domain of parameters 
where the new particles can be found with present day experimental techniques, using upgrades to 
existing experimental facilities. 



Introduction - The standard model (SM) of particle 
physics and theory of general relativity describe correctly 
almost all phenomena observed in nature. Only a handful 
of experimental facts definitely involve particle physics 
beyond the SM: neutrino oscillations (NOs), dark matter 
(DM) , and the baryon asymmetry of the Universe (B AU) , 
which is responsible for today's remnant baryon density 
57b- In Refs. [lj, |2| , it has been suggested that all of them 
may be explained when the matter content of the SM 
is complemented by three right-handed neutrinos with 
masses below the electroweak scale. Different authors 
have investigated aspects of this idea [ll-[3Cij. However, to 
date, it has not been verified that all requirements can 
be fulfilled simultaneously. Claims made in Ref. were 
based on estimates and turn out to be premature from 
today's point of view due to constraints on the properties 
of DM sterile neutrinos from Ly Q and x-ray observations 
that were not known at that time. These constraints 
can be resolved if DM production is enhanced by a lep- 
ton asymmetry generated after sphaleron freeze-out. We 
performed the first quantitative study of this process to 
identify the range of sterile neutrino parameters that al- 
low us to explain at once NO, DM, and the BAU. In this 
letter, we mainly present results; details are given in a 
more detailed publication [3l|- The centrepiece of our 
analysis is the study of neutrino abundances in the early 
Universe from hot big bang initial conditions to temper- 
atures ~ 50 MeV. We combine the results with bounds 
from direct searches for sterile neutrinos and constraints 
from big bang nucleosynthesis (BBN) , which we reexam- 
ined in the face of recent data from neutrino experiments. 
We verify for the first time that right-handed neutrinos 
with experimentally accessible masses and mixings can 
solve all these outstanding problems without any new 
physics above the Fermi scale. We identify the experi- 
mentally interesting parameter region for future searches. 



The Neutrino Minimal Standard Model (vMSM) - The 
scenario outlined above is realized within the ^MSM, de- 



scribed by the Lagrangian 
CvMSM = C-SM + iVR@VR - L L Fi/ R $ - PrF^L L ^' 

-\{v c R M M v R + v R Ml I v c R ). (1) 

We have suppressed flavor and isospin indices. Csm is 
the Lagrangian of the SM. F is a matrix of Yukawa cou- 
plings, and Mm is a Majorana mass term for the right- 
handed neutrinos v R , Ll — (vl, &l) T are the left-handed 
lepton doublets in the SM, and $ is the Higgs doublet. 
We chose a basis where the charged lepton Yukawa cou- 
plings and Mm are diagonal. The Lagrangian JTJ) is well- 
known in the context of the seesaw mechanism [32-35]. 
In the z/MSM, the e igenvalues of Mm are below the elec- 
troweak scale [ll3| . This mass pattern is required to 
simultaneously explain BAU and DM; at the same time, 
it avoids the "hierarchy problem" of the SM in the scale- 
invariant version of the vMSM jH, El] • The z^MSM is 
motivated by the principle of minimality; in comparison 
with the SM, there is no modification of the gauge group, 
the number of fermion families remains unchanged, and 
no new energy scale above the Fermi scale is introduced 

011. 

In the j/MSM, neutrino masses are generated from 
Dirac masses mo = Fv and Majorana masses Mm by 
the seesaw mechanism (v is the Higgs vacuum expec- 
tation value). In the limit Mm 3> tud, there are two 
distinct sets of neutrino mass eigenstates. The block di- 
agonalization of the full mass matrix yields mass matri- 
ces for active and sterile neutrinos m v = —6Mm9 t and 
M N = M M + \(&9M M + Mj^e*), respectively. The 
active mass eigenstates Vi with masses mi are mainly 
mixings of the SM neutrinos vl,o.\ the remaining three 
sterile neutrinos Nj with masses Mi are mainly mixings 
of vrj. Transitions between both are suppressed by the 
active-sterile mixing matrix 9 = mrjM^ 1 . 

In the following, we distinguish between three different 
scenarios. In scenario I, no physics beyond the z/MSM is 
required to explain simultaneously NO, 51dm and J7b, 
as outlined in the Introduction. One sterile neutrino 
(Ni) constitutes all DM. This implies that (a) its mass 
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and mixing are consistent with astrophysical constraints 
and (b) thermal production can account for the observed 
£Idm- The other two (^2,3) produce the BAU and gen- 
erate active neutrino masses via the seesaw mechanism. 
In scenario II, the roles of the Nj are the same, but we 
drop requirement (b), i.e., we assume that DM is made 
of Ni, that was produced by some unspecified mecha- 
nism. In scenario III, we consider the i/MSM as a theory 
of baryogenesis only and drop any constraints related to 
DM. 

To explain the observed DM density J7dm in scenar- 
ios I and II, the lifetime of N% must be larger than the 
age of the Universe. Its decay leaves a distinct x-ray line 
of energy Mi/2 that can be searched for [3^, H(| with x- 
ray satellites. Combining x-ray observations with simula- 
tions of structure formation and phase space arguments, 
it was found that N\ mass and mixing are constrained 
to 1 keV< Mi < 50 keV and 1CT 13 < sin 2 (20 Q i) < 1CT 7 
0. The seesaw relation m, ~ —9Mm9 T implies that the 
coupling of Ni is too small to contribute significantly to 
the active neutrino mass matrix, and one active neutrino 
is effectively massless in the z/MSM 0, [l5[ . 

In scenario I, N\ must be produced thermally in the 
early Universe due to active-sterile mixing [4l|. In the 
absence of lepton asymmetries, (i a = 0, the resulting 
spectrum was found in [23]. For /i Q ^ 0, the Ai dis- 
persion relation in the primordial plasma is modified, 
which results in a resonantly amplified Ni production 
@, 0, E2, S|. This adds a nonthermal, colder compo- 
nent to the Ai momentum distribution, x-ray observa- 
tions, structure formation simulations, and Ly Q forest 
data frl I44T - I46] ] suggest that if \i a = 0, Ai cannot account 
for the observed DM (see figure [lj . Then the presence 
of considerable lepton asymmetries /i Q > 8 • 10 -6 |4j 
becomes a necessary condition for sterile neutrino DM 
production. 

The thermal history of the universe in the z/MSM dif- 
fers from that of the SM, as new interactions generate 
lepton asymmetries /i a ^ during production, oscilla- 
tions, freeze-out, and decay of A/, when all Sakharov 
conditions [52] are fulfilled. No significant lepton asym- 
metries or Ni abundances are created during reheating 
after inflation due to the smallness of F [5j. Baryoge- 
nesis occurs via sterile neutrino oscillations during their 
thermal production 0, E3] in processes as tt — > vA at 
T > T sp h, where the temperature of sphaleron freeze-out 
is T sph ~ 140 GeV for a Higgs mass m H = 126 GeV. Al- 
though the total lepton number violation is suppressed by 
Mi/T <C 1, opposite sign asymmetries are created in the 
sterile and active flavors. The latter are partly converted 
into a BAU by sphaleron processes [54j. To explain the 
observed BAU [55j . a lepton asymmetry [i a ~ 10~ 10 
is required at the sphaleron freeze-out (T ~ T sp h). In 
scenarios I and II, only ^2,3 are produced in significant 
amounts at T > T sp h, as the Ai coupling is constrained 
to be tiny. Soon after, they reach equilibrium and fi a 
are washed out. They exit equilibrium when 11 — » vA 
scatterings freeze out (T ~ few GeV) and decay subse- 
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FIG. 1: Different constraints on JVi mass and mixing in sce- 
nario I. The blue region is excluded by x-ray observations, the 
dark gray region Mi < 1 keV by the Tremaine-Gunn bound 
|47H49l |. For points on the upper solid black line, the observed 
Qdm is produced for /i a — [4]; points on the lower solid black 
line give the correct S7dm for |/i Q | = 1.24 ■ 10 -4 , the maximal 
asymmetry we found at T = 100 MeV. The region between 
these lines is accessible for < < 1.24 • 10~ 4 . Obser- 

vations of the matter distribution in the Universe constrain 
the DM free streaming length. Without resonant production 
= 0), this implies that Mi > 8 keV [44| . which excludes 
the upper black curve and makes resonant production neces- 
sary. Combining both production mechanisms (|^ Q | > 10 ), 
this bound relaxes to Mi > 2 keV [441 ]. However, we do not 
display it here because it depends on ^ Q in a complicated way 
and the calculation currently includes considerable uncertain- 
ties 0, cf. also [13, HI. 



quently (T < 1 GeV). These nonequilibrium processes 
create new lepton asymmetries at a late time. The DM 
production in scenario I is amplified resonantly at tem- 
peratures around T ~ 100 MeV due to the presence of 
these asymmetries. 

We first focus on scenario I. The two requirements (i) 
fi a - 10~ 10 at T ~ T sph (for BAU) and (ii) \fx a \ > 8-10" 6 
at T - 100 MeV (for DM) can be used to constrain the 
properties of sterile neutrinos. Since N% does not con- 
tribute to high and low temperature leptogenesis, this 
can be done in an effective theory with two sterile neu- 
trinos, ^2,3. This theory contains 11 new parameters in 
addition to the SM. They can be chosen as two active 
neutrino masses m, and three mixing angles, a Dirac and 
a Majorana phase, two Majorana masses in Mm and one 
extra complex parameter, associated with CP-violation 
in the sterile sector. 

Condition (ii) implies much stronger constraints than 
(i), so we only consider it in what follows. For the al- 
lowed Yukawa couplings, the asymmetry (ii) can only 
be generated if the CP-violating terms are resonantly 
amplified by a mass degeneracy between M2.3 ~ M 
Q. The asymmetry generation is most efficient for 
r^r ~ H ~ oj, where T^r is the thermal Ni width, H 
is the Hubble rate, and u> is the frequency of oscil- 
lations. It is related to the physical mass splitting SM 
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at the time of low-temperature lepton asymmetry gener- 
ation via uj ~ MSM/T if M < T or simply cj = SM 
if T < M. Due to the interplay of thermal, Dirac 
and Majorana masses, SM is a complicated function of 
M M = diag(M - AM/2, M + AM/2), F and T [3l|. 
The required SM should be smaller than 10 _6 eV for 
M = 2 GeV. Since this is much smaller than active neu- 
trino masses, two unknown parameters of the i/MSM are 
almost fixed by this constraint. For M ~ 2 GeV, the 
splitting of the Majorana masses AM must be equal 
with one part in 10 4 — 10 6 to the mass difference of 
the active neutrinos Q. In addition, the combination 
\Re(rn ! D mD)23\/M must be smaller than active neutrino 
masses by ~ 4 — 6 orders of magnitude [3l| . These tun- 
ings ensure that the Higgs induced contribution to 5M 
cancels the Majorana term AM. Scenario I can only be 
realized within the constrained i/MSM defined by these 
tunings. In the constrained i/MSM, 7 out of 11 param- 
eters are almost fixed either by experimental data or by 
the requirement (ii). The remaining four parameters are 
the common mass M of -^2,3, two CP-violating phases 
in the active neutrino mass matrix, and yet one extra 
CP-violating parameter in the sterile neutrino sector. 



The high degree of tuning SM/M < 10~ 13 in scenario 
I is not understood theoretically. Some speculations can 
be found in (1,0, Ell 56]. However, the origin of this fine 
tuning plays no role for the present work. Scenario II only 
requires the weaker condition (i) that can be achieved 
by a much weaker tuning in AM/M ~ 10 -3 [2l|. In 
scenario III, all three sterile neutrinos can participate 
in baryogenesis. Due to the additional sources of CP 
violation, there is no need for a mass degeneracy (29j . 
Note that this also implies that no degeneracy is needed 
in scenario II if more than three fields ur i are added to 
the SM. 



Method and Results. - The rates of interaction of the 
SM fields exceed those of sterile neutrinos iVj and the rate 
of the Universe expansion. Therefore, the SM sector can 
be described by four numbers: the temperature T and 
the asymmetries p a . The effect of the Nj on the time 
evolution of temperature and on the effective number of 
degrees of freedom is negligible. The state of iVj can be 
described by matrices of density pjv and p^y, commonly 
used in neutrino physics [57j . which allow us to incorpo- 
rate coherences and oscillations between the two flavors 
and are probably sufficient for our purpose, cf. (2j| l57l - 
[76| . The diagonal elements of the 2x2 matrices pn and 
Pjq are the abundances of particles and antiparticles, re- 
spectively, defined as the helicity states of the Majorana 
fields Nj. The time evolution of neutrino abundances is 



governed by the kinetic equations 

1^- = [H,p N }- l -{T N ,p N -p e «}+ l -n a T% , (2) 



dPN 



dT 

■ dp a 

'It 



= [H*,PN\-^*N,PN~p eq }~^ a T a N * ,(3) 



-iYf^pa + itr 



n{ PN -p^) 



-iti 



n*(PN-p eq ) 



(4) 



Here, p eq is the equilibrium density matrix, H is the dis- 
persive part of the finite temperature effective Ni Hamil- 
tonian, and Tpj, T£ and Y a l are rates that are respon- 
sible for dissipative effects, which are calculated from 
thermal field theory [3l|- They describe sterile neutrino 
production, oscillations, freeze-out and decay. Due to 
the various involved time scales, the dependence of the 
asymmetries on model parameters can only be estimated 
under certain assumptions [3], [29|, |6(| [62j that are too 
simplifying for a quantitative study. 

We focussed on scenarios I and II. The most impor- 
tant properties of ^2,3 from an experimental viewpoint 
are their masses M2,?, — M and mixings with active neu- 
trinos. The latter can be parameterized by the quan- 
tity U 2 = tr(#t#). In order to identify the range of M 
and U 2 consistent with conditions (i) and (ii), we calcu- 
lated the lepton asymmetries from T 3> T sp h down to 
T = 100 McV as a function of unknown parameters iden- 
tified above and varied the others within admitted la 
uncertainties. 

Fixing all kno wn neutrino parameters to the values 
given in [77||ll4| . we first identified the Dirac and Ma- 
jorana phases that maximize the produced asymmetries 
at T = T sph and T = 100 MeV in different regions in 
the parameter space. We then scanned for all possible 
values of the remaining parameters. We performed the 
analysis several times with different grids. This allows us 
to identify the parameter regions where condition (i) or 
(ii) or both can be fulfilled. They correspond to sterile 
neutrino properties for which the i^MSM can, along with 
NOs, explain the observed fie , !Jdm or both. We studied 
the mass range 1 MeV < M < 10 GeV; for bigger masses 
it is very unlikely that A^2,3 can be found experimentally 
in the near future. 

The 2/MSM parameter space is also constrained by di- 
rect searches for sterile neutrinos (l7l. I2H l78l - l9Cl | . Here, 
we focus on the most relevant bounds, coming from the 
NuTeV 0, CHARM ^ and CERN PS191 ex- 
periments. The experimental constraints on active-sterile 
mixing have recently been interpreted in the context of 
the seesaw Lagrangian (fT]) [H|,|2(j, cf. also [zl, We 
used bounds on 9 imposed by the negative results of 
[82L 184 , l9lM9a | , provided by the authors of [26J , as input 
and numerically scanned the space of unknown model pa- 
rameters to identify all combinations of M and U 2 com- 
patible with experiment. 

Finally, it is a necessary requirement that ./V2.3 have de- 
cayed sufficiently long before BBN that their decay prod- 
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FIG. 2: Constraints on iVa,3 masses Mi,3 — M and mixing 
U 2 = tr(6'9) in scenarios I (red line) and II (blue line) for 
normal (upper panel) and inverted (lower panel) hierarchy 
of neutrino masses. In the regions within the blue and red 
lines, no physics beyond the z^MSM is needed to explain the 
observed f2a and Qdm, respectively. 

ucts do not affect the abundances of light elements. We 
estimate the inverse iV^ lifetime r by as r -1 ~ itrr^v 
at T = 1 MeV. This is justified as ^2,3 oscillate rapidly 
around the time of BBN. We varied all free parameters 
to identify the region in the M-U 2 plane consistent with 
the condition r < 0.1s. 

Our results are summarized in Figs. [1] and [2l Fig- 
ure [5] shows constraints on ^2,3 mass and mixing coming 
from experiments (green lines) , BBN (dashed black line) , 
neutrino oscillations experiments (black dashed line that 
is labeled seesaw), and cosmology. Scenario II can be 
realized in the region between the blue lines that are la- 
beled BAU. The difference to the result found in [21| is 
mainly due to 813 ^ 0. The region within the red line 
allows us to produce the observed Hdm- It has been de- 
termined for the first time in this work. Although the 
values for the CP-violating phases and AM that maxi- 
mize the efficiency of baryogenesis and DM production 
are very different, the region in which f^B and S^dm can 
be explained simultaneously (scenario I) almost coincides 
with the area inside the red line. In most of the parameter 
space, the relevant CP- violation comes mainly from the 
sterile sector and not from the phases in the Pontecorvo- 



Maki-Nakagawa-Sakata matrix. The masses of A?2,3 are 
correspondingly bounded from below by 1.3 and 2.0 GeV 
for the inverted and normal hierarchies. For the lower 
mass range, sterile neutrinos can be created in decays of 
beauty and charmed mesons, which is crucial for experi- 
mental searches [17J. The two larger eigenvalues of F'F 
can vary between ~ 10~ 12 and ~ 10~ 17 to simultane- 
ously explain Ob and Odm- Typical values correspond 
to Yukawa couplings of ^2,3 of the order 10~ 6 — 10~ 7 , 
smaller than the electron Yukawa by 1-2 orders of mag- 
nitude. The N\ Yukawa couplings are required to be of 
the order 10 -11 — 10~ 12 , smaller than those of by 
5 orders of magnitude, which is comparable to the ratio 
between down and top quark Yukawa couplings. 

Solving Eqs. ©-(lU) allows us to determine the maxi- 
mal lepton asymmetry generated in the j/MSM. Its value 
imposes a lower bound on the mixing of the DM candi- 
date Ni in scenario I. Our result for M < 10 GeV, shown 
in Fig. [1] along with astrophysical bounds, is about 1 or- 
der of magnitude larger than previous estimates [3] . This 
considerably eases the ultimate goal of x-ray searches for 
Nt. 

Conclusion - We performed the first complete system- 
atic study of the j/MSM parameter space, bringing to- 
gether cosmological, astrophysical and experimental con- 
straints. Our results can be summarized as follows: (1) 
Right-handed neutrinos alone can be the common origin 
of neutrino oscillations, DM and the BAU; (2) for a range 
of model parameters, these particles can be found using 
present day experimental and observational techniques. 

The DM candidate iVi can be searched for astrophys- 
ically, using high resolution x-ray spectrometers to look 
for the emission line from its decay in DM dense regions. 
The seesaw partners ^2,3 can either be discovered as 
missing energy in the decay of mesons or by creating 
them in a beam-dump experiment and looking for their 
decay in a nearby detector (l7j . Depending on the 
mass M, different facilities could be used or upgraded 
for this search, including the CERN SPS beam and 
NA62 experiment, LHCb, MINOS, J-PARC or LBNE at 
FNAL. The necessary experiments are challenging due 
to strong constraints on the mixing angle U 2 coming 
from cosmology. 
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